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I. INTRODUCTION 


Presently there is much interest in relativistic 

nucleus-nucleus collisions because of the prediction by quantum 

chromodynamics (QCD) that under conditions of sufficiently high 

energy density a new state of nuclear matter, the quark gluon 

plasma(QGP), can be formed^. The plasma may be observed in real 

time by detecting weakly interacting particles such as photons or 

lepton pairs that emanate from the plasma. Alternatively, it is 

possible to search for a characteristic signal in the angular 

distributions of secondary particles emitted in the collision. The 

latter approach is the only one that can be used when the secondary 

particles are observed in passive detectors such as balloon-borne 
2 

emulsion chambers . 

Recent data from accelerator-based heavy ion experiments 

indicate that compressional nuclear forces are observable in 

3 4 

nucleus-nucleus collisions at high energy ’ . During the collision 
nuclear matter is assumed to behave like a fluid whose 

motion is governed by the equations of hydrodynamics. Some 
fraction of the incident kinetic energy is transformed into 
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potential or congressional energy. The fact that the nuclear 
matter is made up of neutrons and protons is ignored except in the 
final stages of the collision when the matter has expanded and 
cooled. A fraction of the particles emitted in the collision may 
emerge moving in a preferred direction^ allowing for the 
observation of the "footprint" of the QGP production in the angular 
distribution of the secondary particles. 

In this paper we report on the application of statistical 

analyses to the angular distributions of secondary particles 

emitted in relativistic nucleus-nucleus collisions. In particular, 

we have examined a set of JACEE-3 events for the presence of 

correlations in the azimuthal distributions of the secondary 

particles. Two approaches have been used here. In the first case, 

we test for asymmetry in the entire data set of events using the 

7-9 

method of composite unit vectors . Secondly, chisquare 

goodness-of-fit tests and discrete Fourier transform analysis are 
used to examine event-by-event for correlations between particle 
intensity and azimuthal angle. Our analysis is further extended to 
include searches for evidences of non-statistical fluctuations in 
the pseudorapidity distributions of single events. It is expected 
that the investigation of non-statistical fluctuations in single 
events will be most productive for high multiplicity collisions 
since it is assumed that nucleus-nucleus collisions are made up of 
many incoherent elementary processes^ - ^. 
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II. METHOD OF COMPOSITE UNIT VECTORS 


The method of composite unit vectors permits one to test for 
azimuthal asymmetry in the distribution of an entire data set of 
events. A summary of the reactions included in this set and the 
range of energies covered is provided in Table I. In the method, 
the direction of the emitted particle is projected onto the 

A 

azimuthal plane and a unit vector, r^, is constructed from this 
projection. As shown in figure 1., R i describes the emitted 
particle’s direction and R^SIN0^ is the magnitude of its projection 
onto the azimuthal plane. The unit vector in the azimuthal plane 
simply becomes 


r. = (COS0.) i + (SIN0.) j . (1) 

For each event we construct a quantity, w, given by 


w = 



i - l.n , (2) 


n 

where n is the number of secondaries and 0^ is the azimuthal angle 
associated with the ith emitted particle. We note that events with 
a small degree of asymmetry will have w values near or equal to 
zero and that w will increase with increasing asymmetry. If it is 
assumed that the emitted particles are randomly distributed in the 
angle 0, the situation is mathematically equivalent to the problem 
of calculating the probability distribution associated with n 
frequencies of unit amplitude each possessing a random phase. The 
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Table 

I. 

Summary of 
energy for 
available 

reactions treated 
at least one event 
at the time of this 

in this work, 
in this group 
analysis . 

< a) The 
was not 

REACTION 

RANGE OF 
MULTIPLICITIES 


ENERGY RANGE 
(GeV/nucleon) 

N i 
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28 
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Pb 
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Figure 1. Coordinate system diagram 
unit vector method. 

illustrating the composite 



resulting probability distribution, when expressed in terms of the 
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parameter w, is approximated by the exponential form , 

P(w)dw - exp(-w)dw. (3) 

It therefore should be possible to obtain a measure of the 
azimuthal asymmetry in the entire data set by constructing the 
distribution of events in w and comparing its shape to the 
exponential. Deviation from this exponential behavior may then be 
interpreted as an indication of anisotropy in the data. 

The distribution in w for the JACEE-3 events analyzed is shown 
in figure 2. The solid curve in the plot is the renormalized 
exponential function. The dashed line is the result of of a least 
squares fit to a single exponential, 

n = Aexp(-Bw) , (4) 

where A = 24.5694 and B = 0.84962. In general, the agreement 

between the experimental distribution and the theoretical curves is 
good. This result is not unusual in that purely statistical 
fluctuations that occur event to event will be lost when averaging 
over a sufficiently large set events. It is therefore important 
that single events be investigated for evidences of non-statistical 
fluctuation. In the next section of this report we present a method 
for the analysis of structure observed in the azimuthal 
distributions of single events. 
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EVENTS PER BIN 



Figure 2. Asymmetry distribution plot for the 63 JACEE-3 events 



III. SINGLE EVENT ANALYSIS 


A. Chisquare Tests 

In order to examine more closely the question of nonuniformity 
in the azimuthal angular distributions, each distribution was 
subjected to a null hypothesis test wherein the probability that 
its intensity pattern could be attributed to chance occurrence was 
determined. In effect, a least squares fit of the data to the 
average of the distribution was performed and a chisquare value 
given by 



1 A 2 

was evaluated for each event 1 . In the present case X corresponds 
to the sum of the squares of the residuals with each residual 
weighted by its fractional standard deviation. The integral 
chisquare probability^ , Py_ , was also determined for each event. 
It represents the probability that a randomly chosen set of N data 
points will yield a value of X as large as or larger when compared 
with the parent distribution. An assumption that is inherent to 
this part of the analysis is that the experimental X values are 
characterized by a chisquare distribution of degrees of freedom. 
This will indeed be the case if the X^'s are independent and 
identically distributed variables each following a normal 
distribution with a mean of zero and unit variance. The exercise 
of evaluating this assumption is complicated by the relatively 


8 



small number of events in the sample and the fact that the events 
generally correspond to different reactions rather than a collection 
of events from the same reaction. Further, variations in the 

multiplicity among events lead to variations in the values 
associated with the individual angular distributions. Thus, the 
events that fall within a given bin of the experimental X 2 
distribution will in general have a range of values. With these 
considerations in mind we have taken the following approach to 
obtaining a comparison between theory and experiment. 

The general form of the chisquare distribution is given by 




a¥ y - 2W eM-7 2 n) 


( 2 ) 


T7T^ 


V/2) 


( 6 ) 


where v is the number of degrees of freedom associated with the 

2 

distribution. For each bin in the experimental X distribution we 
construct an average probability given by 


P 


k 



(X-) ( 2)/2 exp(-Xj/2) 

(2) i/2 f( K/2) 


(7) 


where 2 

N = number of X^ values that fall within the kth bin, 

k = index specifying the bin, 

2 

X • = chisquare value obtained for the ith angular 
1 distribution 

and 

y . = number of degrees of freedom for the ith angular 
1 distribution. 
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The values are normalized to satisfy the criterion, 


N_ 


k - 1 1 Pt = "* 


( 8 ) 


2 

is the number of bins in the experimental distribution and 
Ng is the total number of events in the data set. The comparison 
between the above theory and the experiment is provided in figure 
3. The shape of the two distributions is seen to be in good 
qualitative agreement thereby supporting the assumption upon which 
the integral probabilities are based. 

B. Discrete Fourier Transform Analysis of Azimuthal Distributions 

Three events with relatively high statistics were chosen for 
discrete Fourier transform (DFT) analysis 16,17 . The events are 
identified in Table II, and their distributions are shown in 
figures 4a-4c and 5a-5c. It should be noted that the number of 

bins in set 4 is based on an optimization of the sensitivity of 
the null hypothesis test. The criterion for the Fourier analysis 
is that the number of bins be an integral power of two. The exact 
discrete transform was determined for each distribution according 
to the relationship 

V 2 

N 1 ■ N av + l \[COS(k„t +/0 )] , (9) 

k = 1 


where 



k 

A k 

*k 


■ number of bins in the distribution, 

= average number of particles per bin, 

= index over which the sum is performed, 

= calculated Fourier amplitudes, 

= calculated Fourier phases, and 

= predicted value of the number of counts in the bin. 
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EVENTS PER BIN 



PARTICLES PER BIN 



PHI (DEGREES) 


Figure 4a. Azimuthal distribution of emitted particles for event FE70642 
Ti — *-AgBr at 41 GeV/nucleon. (1 - P ) = 0.978406. 


PARTICLES PER BIN 



Figure 4b. Azimuthal distribution of emitted particles for event FE6869B, 
Fe— ►O at 55 GeV/nucleon. (1 - P „) = 0.999953. 
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Figure 4c. Azimuthal distribution of emitted particles for event FE61927 
Fe — ►AgBr at 24 GeV/nucleon. (1 - P ) = 0.29763. 
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PHI (DEGREES) 


Figure 5a. Azimuthal distribution for event FE70642. The solid curve 
drawn through the points is the principal component of the 
fourier transform and corresponds to the function 
N i = N a» + A 6 C0S(6 u. tl ♦ d 6 ). 
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PARTICLES PER BIN 



Figure 5b. 


Azimuthal distribution for event FE6869B. The solid curve 
drawn through the points corresponds to the function 
N. = N av + A 1 C0S( cot. + 0 X ) + A 2 C0S(2w t- + 0 7 ). 


PARTICLES PER BIN 



Figure 5c. Azimuthal distribution forevent FE61927. The solid curve drawn 
through the points corresponds to the function 


N, = 


av 


A^COS ( wt.j 


+ 0j) 
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Table II. Reactions for which discrete Fourier transform analysis 
was performed on the azimuthal angular distribution. 


EVENT 

REACTION 

ENERGY 

(GeV/nucleon) 

N 

(1-P ) 
exp* 

FE70642 

Ti > AgBr 

41 

283 

0.978406 

FE6869B 

Fe > 0 

55 

132 

0.9999 53 

FE61927 

Fe > AgBr 

24 

165 

0.297630 


The quantity ojt^ is the angle at which the are calculated. The 

complete results of the Fourier analysis are given in Appendix 

A. Although eight Fourier amplitudes were calculated by the 

procedure, each is not necessarily of equal importance to the 

reconstruction of the original distribution. In order to judge 

the relative importance of a given amplitude, an integral 

chisquare probability, P , has been calculated for each. 

X 

The P is given by 


2 


P 

X 


(xj; v) = f 

O 


*k 


( v ¥ v - 2)/2 



exp(-y 2 /2)dy 2 # 
r(v/2) 


( 10 ) 


where 



.< b i - N ik> 



(U) 


and k refers to a single amplitude in equation (9). According 
to equation (10) large values of P indicate a poor fit 
and small values indicate a good fit* This approach per- 
mits one to easily determine whether the character of a given 
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distribution is dominated by a subset of amplitudes. With this 

idea in mind, the solid curves shown in figures 5a-5c were drawn to 

show those amplitudes that have the smallest probabilities of 

chance occurrence. Of the three events, FE6869B appears to be the 

18 

only one well characterized by two slowly varying components . The 
distribution for event FE70642 on the other hand is best 
characterized by a rapidly varying amplitude which is more typical 
of noise, 

C. Fluctuation Analysis of the Pseudorapidity Distributions 

19 

We have used the method of Takagi to test for nonstatistical 

fluctuations in the pseudorapidity distributions of individual 

events. This method involves several steps that are only 

summarized here. First, a smooth function, f s (7J) that reflects the 

average properties of the distribution is determined by least 

squares fit of the experimental distribution to a superposition of 

20 

Chebyshef polynomials . The number of terms in the sum is 
governed by the expansion that minimizes the chisquare value for 
the fit. Thus, we make the definition 

m 




a 0 /2 



a k T k (z) =f^(7)). 


( 12 ) 


If there are L bins in the histogram of the distribution, 

m<L-l and z = ~7] max )/( T) max ' ^ min h T k (z) 

is the Chebyshef polynomial of degree k and f('Tj) is 
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known for the range Tj 
V, of the experimental histogram is then used as a measure of the 
fluctuations in the distribution, i.e., 

L+l 

V exp ' Y, If<,11l) ' f( ' 7?i - l) ' ll3> 

i=l 

where f(7^ Q ) = f (T^ + ^)=0. The pseudorapidity distribution is then 

simulated by a Monte Carlo procedure. In the present case the 

21 

method of Von Neumann rejection was used to produce the 
pseudorapidity distributions for many "equivalent” events. For 
each simulation, the number of secondary particles was held fixed 
to the experimentally observed value. Points within the range 
^max— min were sam pled according to f g (7 ~j). The Monte Carlo 

data obtained should contain distributions of purely statistical 
fluctuations that can be compared with the fluctuations in the real 
events. 

Two examples of smooth fits, f g (7 ]), are provided in figures 6 

and 7 together with their empirical distributions. For comparison, 

the distributions resulting from single event simulations are shown 

in figures 8a and 8b for events FE6869B and FE70642, respectively. 

As expected, distributions based on the average of 250 event 

simulations reflect the features of the empirical fits for these 

events. The plots are provided in figures 9a and 9b. Note further 

from figures 10 and 11 that the experimental fluctuation values, 

V , fall well within the simulated V distributions thus 

exp 

indicating that the fluctuations in the empirical data are largely, 
if not purely, statistical in nature. 


max —'V— ^min* The total vertical length. 


20 





dN/dn 


240 


220 


FE 70642 

Ti — ►AgBr at 41 GeV/nucleon 



Figure 7. Experimental pseudorapidity distribution for event FE70642. 

The solid curve is the result of a least squares fit of the 
data to the function: 

0 

f s <"> = v 2 + y\ T k ‘ z >- 

k=l 
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Monte Carlo Data 

Ti — ►AgBr at 41 GeV/nucleon 

1 event 


N s * 283 



Figure 8b. Single event simulation of the distribution function, 
fg(n), for event FE70642. 
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Monte Carlo Data 

Fe — at 55 Gev/nucleon 

250 events 


N s - 132 


u 



2,1 


T 

v 



Figure 9a. This plot shows the average distribution resulting from a 

4 

250 event simulation of the function, f (n)» for event 
FE6869B. 
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Figure 9b. This plot shows the average distribution resulting from a 
250 event simulation of the function, f® ( n ) , for event 
FE70642. 



Number of Events 



Figure 10. Distribution of the theoretical fluctuations 
V- for event FE6869B. 




of Events 



Figure 11. Distribution of the theoretical fluctuations 
V^., for event FE70642. 


IV. WOUNDED NUCLEON MODEL SIMULATIONS 
OF NUCLEUS-NUCLEUS COLLISIONS 


22 

In the wounded nucleon model it is assumed that the 

inelastic collision of two nuclei is an incoherent composition of 

collisions of individual nucleons and those nucleons that undergo 

at least one inelastic collision are termed "wounded". The size of 

the fragment and the number of alpha particles that are emitted in 

the collision was determined by random sampling of empirical data 
23 

tables . The remaining mass involved was then divided into wounded 
and unwounded nucleons. The wounded fraction was found by sampling 
a distribution which was a binomial approximation with an average 
value of 0.4. We used this model to simulate events FE6869B and 
FE70642 and test the azimuthal angular distributions to see how 
many possessed (1 - P^) values as large as or larger than the 

experimental values. Both events were simulated 100 times. In each 
case a relatively small fraction, -13%, of the simulated angular 
distributions had (1 - P^) values satisfying this condition. 

V. SUMMARY AND CONCLUSIONS 


The experimental investigation of heavy nuclei collisions at 
high energies is necessary for establishing to what extent 
collective phenomena are important in nuclear collisions and to 
identify those aspects of elementary collisions that get amplified 
in the nuclear interactions. The methods presented in this paper 
enable one to place on a quantitative basis the decision as to 
whether distributions of particles emitted in a collision exhibit 
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nonrandom structure and to characterize the structure in terms 


of its periodic behavior. In particular, we have seen that 
discrete Fourier transform analysis is quite effective at the 
latter. We also find that the present data set supports the 
assumptions upon which our statistical analyses are based. 

One event, Fe >0 at 55 GeV/nucleon shows an azimuthal 

distribution with a bimodal character describable in terms of a 
two-component Fourier representation. Results of computer 
simulations, based on the wounded nucleon picture of the collision, 
are consistent with the large deviation from uniformity observed in 
the distribution for this event. However, we do not observe 
evidences of non-statistical fluctuations in the pseudorapidity 
distribution for this or other events in the data set. 

In order to take on physical meaning the results of these 
statistical analyses must be combined with deterministic modelling 
based upon the mechanisms in effect during collision. In this 
regard, it is important to extend the present activity to include a 
hydrodynamical description of the collision that takes account of a 
high energy density and small nuclear radius. This work is 
presently in progress. 
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